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[57] ABSTRACT 

A deconvolution approach to adaptive signal processing has 
been applied to the elimination of signal multipath errors as 
embodied in one preferred embodiment in a global position- 
ing system receiver. The method and receiver of the present 
invention estimates then compensates for multipath effects 
in a comprehensive manner. Application of deconvolution, 
along with other adaptive identification and estimation 
techniques, results in completely novel GPS receiver archi- 
tecture. 
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METHOD AND APPARATUS FOR 
REDUCING MULTIPATH SIGNAL ERROR 
USING DECONVOLUTION 

ORIGIN OF THE INVENTION 

The invention described herein was made in the perfor- 
mance of work under a NASA contract, and is subject to the 
provisions of Public Law 96-517 (35 U.S.C. § 202) in which 
the Contractor has elected not to retain title. 

TECHNICAL FIELD 

The present invention relates to signal error reduction, 
and more specifically to a method and apparatus for reduc- 
ing Global Positioning Satellites (GPS) multipath signal 
error using deconvolution to determine the attitude and 
position of kinetic platforms. 

BACKGROUND ART 

The superiority of GPS receiver technology over other 
competing technologies, in terms of the cost, reliability, 
mass, size, and power considerations, as proven in its 
various applications such as navigation, spacecraft orbit 
determination, and surveying, has constantly prompted the 
extension of GPS to other important areas that require 
increasingly more exacting performance from the GPS sys- 
tem. One such area, for example, is the application of GPS 
to the attitude determination of aircraft and earth orbiting 
satellites. 

For the GPS receivers to provide the requisite precision, 
the existing error sources in the GPS system must be 
eliminated or greatly reduced. At the current state of GPS 
technology, the most significant error source is the signal 
multipath propagation. For example, in spacecraft opera- 
tional environments, the GPS signal is reflected from vari- 
ous structural components of the spacecraft and these 
reflected components are received by the GPS receiver along 
with the desired direct line-of-sight (LOS) path. The 
reflected signals differ from the desired LOS path signal in 
terms of their delays, amplitudes, and phases. The carrier 
phase tracking loop provides no inherent discrimination 
against the multipath signals and thus tracks the phase of the 
composite signal corrupted by multipath components. The 
resulting differential carrier phase estimation error can be 
orders of magnitude higher compared to the case of no 
multipath propagation in many GPS applications. For 
example,, the measurements obtained by the RAD CAL 
satellite GPS-ADS (Attitude Determination System) experi- 
ments have shown that the differential range error in such 
environment is of the order of 1 cm corresponding to an 
attitude determination error of about 0.5 degree. Thus for the 
GPS receivers to provide precision pointing knowledge 
(order of 1 arcmin or better with 1 meter antenna baseline) 
or a differential range accuracy of about 0.3 mm or better, 
the multipath effects must be suppressed by orders of 
magnitudes. Similar accuracy may also be desirable in other 
GPS precision applications in the presence of multipath 
signals such as GPS based geophysical measurements and 
precision surveying. 

Among the past approaches to deal with the multipath 
problem, one approach involves reducing the early -late 
delay spacing among the correlators in the GPS receiver 
code lock loop. However while this approach reduces the 
code range errors to some extent, it does not aid in the carrier 
phase measurements accuracy that is the basis of most GPS 
precision applications. Moreover even the reduction in the 


2 

code range error is limited and if the early -late spacing is 
smaller than the initial delay error due to multipath signals 
(easily the case with many multipath situations), then the 
loop error can be very high and the loop may not even track. 
5 In the work of VanNee et al., a set of nonlinear implicit 
equations are derived on the basis of maximum likelihood 
estimation theory. These equations contain the amplitude,, 
phase and delay of multipath signals. The paper proposes to 
solve these highly nonlinear equations in a recursive 
10 manner, but there is no explicit solution presented. The 
method has been implemented in a Nov Atel receiver using 
multiple correlators. Although this approach is useful for 
code phase measurement and a benign multipath 
environment, it may not be suitable for severe multipath 
15 environment and for precision GPS applications based on 
carrier phase measurements. 

STATEMENT OF THE INVENTION 

The present invention provides a method and apparatus 
20 for simultaneous estimation and compensation of the mul- 
tipath estimation errors in both the carrier phase lock and the 
code delay lock loops in GPS receivers. The preceding is 
accomplished by disclosing new techniques for dealing with 
the multipath problem and includes these techniques in an 
25 improved GPS receiver. These techniques are based on the 
application of the optimal deconvolution approach in a 
somewhat unconventional manner, as compared to its appli- 
cation in other fields such as seismology and telecommuni- 
cations. The proposed method consists of first estimating the 
30 impulse response of the effective multipath channel by a 
least squares algorithm. This step is followed by obtaining 
an inverse filter which equalizes the multipath channel 
response to the desired ideal multipath free response to the 
maximum extent possible within the specified constraints of 
35 the implementation complexity. From the equalized 
response, estimates of the true carrier phase and code delay 
can be made. The simulation results demonstrate that the 
proposed method is capable of reducing the multipath dis- 
tortion by several orders of magnitude in an environment 
40 that is more severe than may possibly take place in any 
realistic precision GPS applications environment. 

The invention also is embodied in a receiver comprising 
a series of correlators which generate demodulated signals 
from known reference signals, and the demodulated refer- 
45 ence signals are then used to develop a discriminator func- 
tion from which multipath error can be approximately 
determined. The present invention then compensates for the 
multipath error, and produces a final differential phase and 
code range estimate after compensating for the multipath 
50 error. 

The present invention presents an optimum solution to the 
inherent tradeoffs among hardware/software implementation 
complexity, the extent of the multipath expected in the 
application, and the degree of multipath cancellation. 

55 BRIEF DESCRIPTION OF THE DRAWINGS 

The objects and features of the present invention, which 
are believed to be novel, are set forth with particularity in the 
appended claims. The present invention, both as to its 
60 organization and manner of operation, together with further 
objects and advantages, may best be understood by reference 
to the following description, taken in connection with the 
accompanying drawings. 

FIG. 1 is a block diagram of the elements of a receiver in 
65 a preferred embodiment of the present invention; 

FIG. 2 is a flow diagram of the channel impulse response 
estimator and equalizer of FIG. 1; 
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FIG. 3 is a graph illustrating the tracking of the Discrimi- 
nator Function amplitude vs. delay in simulation example 1; 

FIG. 4 is a graph illustrating the Discriminator response 
vs. delay in simulation example 1; 

FIG. 5 is a graph of the normalized delay error magnitude 
vs. iterations for simulation example 1; 

FIG. 6 is a graph of the pointing error vs. iterations for 
simulation example 1; 

FIG. 7 is a graph of the pointing error vs. iterations for 
simulation example 1 with Delay >0.9; 

FIG. 8 is a graph of the normalized delay error magnitude 
vs. iterations for simulation example 2; 

FIG. 9 is a graph of the pointing error vs. iterations for 15 
simulation example 2; 

FIG. 10 is a graph of the pointing error vs. iterations for 
simulation example 3; and 

FIG. 11 is a graph of the normalized delay error vs. 20 
iteration number for simulation example 3. 


4 

In the following, the deconvolution approach is extended 
to the problem of multipath elimination in the GPS receiver 
code tracking and carrier phase lock loops. Referring to the 
drawing of GPS receiver multipath compensator 100 shown 
5 in FIG. 1, the antenna 1 receives the GPS signal which is 
filtered and amplified by the RF front end 2. In the absence 
of the multipath signals, the output 3 of the RF front end 2 
is given by: 

= A c cos[w c r + a(r)-J ® 

where the receiver noise is not considered in the first 
instance, A c is the received signal amplitude, oo c is the GPS 
signal frequency, and a(t)=±l is the pseudo-random code 
waveform that phase modulates the carrier. It is assumed that 
the data modulation is removed in a decision-directed man- 
ner. In the presence of N multipaths in addition to the direct 
line-of-sight path, the input signal may be characterized as 

s m (t)=4t)+ a At-<)+- ■ ■ +a jV s(f-T iV ) (3) 


DETAILED DESCRIPTION OF THE 
INVENTION 

The following description is provided to enable any 
person skilled in the art to make and use the invention and 
sets forth the best modes contemplated by the inventors of 
carrying out their invention. Various modifications, 
however, will remain readily apparent to those skilled in the 
art, since the generic principles of the present invention have 
been defined herein specifically to provide a method and 
apparatus for reducing multipath error especially in GPS 
applications. 

In the more conventional telecommunication applications 
of the deconvolution approach or the equalizer theory, the 
multipath propagation channel is modeled as, 


where a,- and x i denote respectively the amplitude and delay 
of the i th multipath for i=l,2, . . . ,N, and s(t) is the direct 
path or line-of-sight component of s m (t). After substituting 
(2) in (3), the composite received signal may be expressed 
as follows: 

s m (t) = A c cos^co c t + a(r)-) + A c cos(w c r + -a(t-T l ) + 9 m{ j + 

...+a N A c cos(w c r + ^a(t - t n ) + 0 mN ); 9 m . = -w c r ; 

This signal s m (t) is translated from the carrier frequency 
to an intermediate frequency (IF) by converter 4 to yield IF 
signal s /i7 (t) at an intermediate frequency coj. Thus, equation 
(4) can then be represented as 


25 


30 


35 


« (1) 
yj= h k u j-k + vj 40 

*=-<71 

where {uy} represents the transmitted symbol sequence, {y 7 } 
is the channel output sequence, and {h , h +1 , . . . , h 0 ,h l7 45 
. . . h q2 } represents the discrete channel impulse response. 
The additive noise sequence {vy} is usually assumed to be a 
zero-mean white Gaussian. The basic process of deconvo- 
lution involves the estimation of the transmitted input 
sequence {%} on the basis of noisy observations {y^} 50 
assuming that the discrete channel impulse response {h^} is 
known to the receiver. In the case of unknown channel 
response, adaptive equalization techniques may be used 
wherein first an approximate estimate of {h^} is obtained on 
the basis of a training sequence known in advance to the 55 
receiver and the channel output {y*} and subsequently the 
estimate of {h*} is refined adaptively with u k replaced by its 
estimated/detected version in the adaptive algorithm. From 
the real time estimate of {h*}, a time-varying inverse filter 
is derived which then filters {y*} to obtain the estimate of 60 
{%}. In practice the two steps of estimating {h^} and then 
finding the corresponding inverse filter are combined into a 
single step of finding directly the adaptive equalizer coef- 
ficients. The problem wherein the adaptive equalization is 
achieved without any training sequence is relatively more 65 
difficult and has also received considerable attention in the 
literature. 


s m (t) = A c cos ( w 1 f + a ( r)-j + ar 1 A c cos ( c t > 1 r + - a ( r - r 1 ) + 0 J „ 1 ) + 

...+a N A c cos ( w ! t+^a(t-T N ) + 9 mN ); 9 m . = - a) c T h 

In the conventional delay -lock discriminator the signal in (5) 
is correlated with the reference signals 

/ 7 T \ (g) 

s L (t) = 2 cos ( w ! t + a(t - t - T d ) - + 9 C J 
s E (t ) = 2 cos ( o >! t + a(t - T + T d )~ + 

where t is the delay tracking error, x d is the offset delay 
referred to as early and late correlator delay and 0 C is the 
reference oscillator phase. The resulting correlation func- 
tions are given by 

% (9^L(9= fl ca(9-a(^ -^d) COS ( 0 /»n- e c)+- ■ ■ +KV-N 

a(f-%).a(r-T-Td)-cos(6^-6 c ). (7) 

The resulting signal denoted by R ll (t) may be written in the 
form below. 

R^)=Kfic^d)+h 1 ^Jx+x d -x 1 )+. . . +h Ni R c (x+x er x N )h k ~ 

a * 4 c cos ( 0 rajt - 0 c ); fe = O , l , . . . ( 8 ) 

In equation (8) R c (t) represents the code autocorrelation 
function and a 0 =l. Similarly the corresponding signal 
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obtained with s L replaced by S E in (7) and denoted by R 2l (t) 
may be written in the following form. 


R 2i {x)=h 0i R c (T-T d )+h li R c (T-T ( rT 1 )+. . . +h N fi c (x-x d -x N ) (9) 

The difference between R lt - and R 2i , the so called discrimi- 
nator function D/x) is now given by 

D i (r)=h 0igc (x)+h ligc (x-x 1 )+. . . +h Nigc (x-x N ) (10) 

where 

^ c (T)=i? c (T+T)-i? c (T-T rf 

is the discriminator function in the ideal case. Clearly if 
a A =0 for all k^O then D,-(x)=g c (x) corresponding to the case 
of no multipath propagation. The standard delay lock loop 
converges to the solution x m of the equation D/x)=0 instead 
of converging to 0 which is the solution of g c (x)=0. The 
multipath error x m can be excessive depending upon the 
actual multipath environment encountered as will be illus- 
trated by several simulation examples. 

In the approach of the present invention the discriminator 
function is measured first and then with equation (10) is used 
to estimate all the unknown variables including the multi- 
path delays, amplitudes, and phases and the estimates of the 
multipath errors both in the phase lock and code lock loops 
via deconvolution. The multipath errors are then compen- 
sated for in arriving at the final differential phase and code 
range estimates. 

To accomplish these objectives, multiple correlators cor- 
responding to multiple values of x d are used in the present 
invention instead of just two correlators. Returning to FIG. 
1, signal 5 is correlated with carrier signal 43 Q(t) in 
correlator 6, and the output 7 is correlated in correlators 
8,12,16 with signals 54, 53, and 55, respectively. The 
outputs 9,13, and 17 are received by integrate/dump units 
10, 14, and 18 and each output is integrated over T seconds 
to generate the correlator functions represented in equation 
(8). In the same manner, all the correlator functions R^x+iA) 
for all integers i in the range of -Q 1 to Q 2 are generated using 
additional correlators and integrate/dump units. Output 11, 
R,(x+Q 2 A) is the same as R n (x) of equation (8) with x d in 
(8) replaced by Q 2 and the carrier signal 43 C,-(t) is given by 

Ci(t)= 2cos ((Ojf+Oj. (11) 

The multipath resolution parameter A is selected sufficiently 
small to provide the required multipath resolution and the 
number of correlators ((Jj+Q 2 +1) is dependent on the extent 
of multipath delay spread and A. The linear combiner 20 
receives output 11, 15, and 19 from the respective integrate/ 
dump units 10,14, and 18, as well as any other integrate/ 
dump units (not shown) and generates the sampled values of 
the function D,-(x) given in (10) by differencing pairs of its 
input, that is, 

D i (x+kA)=R i (x+(k+I d )A)-R i (x+(k-I d )A);k=-M i , . . . , 

0 ,M 2 . (12) 

where x d is selected equal to I^A for some integer l d and M 1 
and M 2 values are dictated by the extent of multipath. Thus, 
in terms of M 1 and M 2 , the Q 1 and Q 2 are given by 
Q 1 =M 1 +I^, Q 2 =M 2 +I^. The output 36 from the linear com- 
biner 20 D; is the (M a +M 2 +1) vector whose k th component 
is given by D/x+kA) in (12). 


6 

The signal 5, S /i7 (t) is correlated in 21 with the carrier 
signal C^(t) 45 which is the quadrature phase version of C/t) 
and is given by: 

5 

C q {t)=- 1 sin ((Ojf+Gj (13) 

Correlating output 22 of the correlator 21 with signal a(t-x) 
in correlator 27 and integrating the output 28 over an interval 
of T seconds in the integrate/dump unit 29 generates R^x) 
10 which is related to the parameters of interest by equation 
(14) as 

■ ■ +/^cCt-Xv) /^= a*A c 

sin(e mr 0 c );k=O,l, . . . # (14) 

In the same manner all of the correlator functions R^x+iA) 
for all integers i in the range of -Q 1 to Q 2 are generated using 
additional correlators exemplified by 23,31 and whose out- 
puts 24,32 are received by integrate/dump units 25,33 as 
2Q shown in FIG. 1. Outputs 26,30,34 are transferred to the 
linear combiner 35 which generates (M 1 +M 2 +l) vector T) q 
whose k th component D^(x+kA) is given by: 

D q (x+kK)=R q (x+(k+I d )A)-R q (x+(k-I d )A; k=-M 1} . . . , 

0 ,M 2 . (15) 

25 

From equation (15) it follows that D^(x) may be expressed 
as: 

D q (r)=h 0qgc (x)+h lqgc (x-x^+. . . +h Nqgc (x-x N ). (16) 

30 For the purpose of estimating and filtering, the outputs 36 
and 37 are combined into the vector D given by 

D=D i+ jD q (17) 

35 

where j= v -l and the k th component of D denoted by D* is 
given by 

D k =D(x+kA)=D i (x+kA)+jD q (x+kA) (18) 

40 From equations (10), (16), and (18) it follows that 

D(x)=h ogc (x)+h lgc (x-x 1 )+. . . +h Ngc (x-x N )+n(x) h^h^+jh^; 

k=l, . .. ,N. (19) 

Here, D(xk) is the complex-valued discriminator function, 
45 and n(t) represents the noise at the correlator output corre- 
sponding to the noise at the receiver input. 

FIG. 2 illustrates a flow chart of the steps executed by the 
channel impulse response estimator and equalizer 38. Step 
200 involves obtaining the two discriminator functions D, 
50 and T> q from the linear combiners 20,35 and solves for 
possible channel impulse responses h* based on delay kA. 
The output from step 200 is a series of channel impulse 
responses h~ N through m, which are then processed in step 
210. Step 210 selects the impulse response which best fits 
55 the measurement for the inputs D, and D q , denoted as m°. In 
step 220, the corresponding deconvolution filter coefficients 
j are computed from the results of step 210, and the 
coefficients are used in step 230 to remove the multipath 
contribution by convolving measurements {D 7 } with {jy}. 
60 The resultant discriminator function D eq is then solved in 
step 240 for Re{D e? (x)}=0 by linear interpolation of the 
sampled version of D eq . The delay obtained, x is used in step 
250 to solve for the channel impulse response h f and to 
compute the argument (phase) of the central coefficient h 0 f 
65 in the final step 260. 

Turning to the preferred solution, in FIG. 1 estimator 38 
processes the input vectors D, and to estimate the 
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multipath channel impulse response, obtain the coefficients 
of a deconvolution filter, and from these derive an estimate 
of the delay t and the phase difference 0=(9 mO -0c) such that 
these estimates are not corrupted by the multipath signals. 
For an understanding of the estimator 38, Equation (19) is 5 
rewritten in its discrete form (20), that is, substitute x=jA 
where j is an integer into equation (19): 


. . . ,z 0 , . . . ,z (q2:t . K2) }. Therefore the equalized 
discriminator response obtained by convolving {yy} with 
{jy} sequence is given by 

(<72+^2) (23) 

D eq (T)= Yj Zigc(r - /A) 

*'=-(<71 +*1 ) 


<72 _ (20) 

yj= X h kgj- k +nj; j = 

*=-<71 

Here y,=D(jA); g,=&(jA); n~n(jA); and (-M 1; M 2 ) repre- 
sents the interval over which the measured discriminator 
function is significant. Note that equation (20) is of some- 
what more general form than (19) in that this form may also 
include negative values of the multipath delays. The delays 
t-l, . . . , Tjy are approximated by integer multiples of A which 
can be achieved by selecting A sufficiently small. A least 
squares estimation technique is now applied to estimate the 
channel impulse response coefficients {h*} and the corre- 
sponding deconvolution or equalizer filter coefficients from 
{y,}- Letting 


h r =[h_ gi , ■ ■ ■ , h_ 1; h 0 , h ls . . . , h^J 

Sj = \.%j+qi> ■ ■ ■ j 8y+iJ 1> ■ ■ ■ » j =— • • • , 

M 2 

the least squares estimate of the unknown channel impulse 
response vector is given by 


m 2 V 1 m 2 

Y 8jgj+oil Y Wi 

J=-My ) j=~M l 


( 21 ) 


where cr„ 2 =E[n- 2 ] and I is the identity matrix. Now denoting 
by % 0 the truncated (or zero padded depending upon whether 
qp>K; or vice versa; i=l,2) version of h, i.e., with 

Xo = [hjq, • • • A, • • • h-A: 2 ] 

where and 1^ are the selected numbers of the deconvo- 
lution filter coefficients, and by Xj the j times shifted version 
of Xo, 

X/ = [h;q+/'> ... ,hj, , h _ Kl+ j\ 
the optimum parameter vector f is given by 

( (^2 +<72) fj-2 V 1 ( 22 ) 

/= Y xj*j + -7' 

(j=-{K l+qi ) 


Note that in the definition of Xj the entries for hy are set equal 
to zero if j falls outside the interval [-qi,q 2 ]- Also the 
constant k is given by 

J L 2 

/f= — V g)\ l = l 1 +l 2 + i 

j=-l i 


4- q i ; L,2 = M 2 + <%2' 


The equalized channel response z is obtained by convolving 
the two sequences {h_ qi , . . . ,ho, . . . ,h ? J and ■ ■ ■ 

,j 0 , ■ ■ ■ Jk 2 \- The elements of z are denoted by 


For the case of perfect equalization z 0 =l, z,=0 for i^O and 
D e? (T)=g c (x). In practice Re{D e? (x)}«g c (x) and the solution 
of the equation 

Re{D eq ( x)}-0 (24) 

15 is the estimate of the delay x. When the initial delay 
uncertainty x is smaller compared to A/2 in magnitude, the 
j th components of the vector D, and T) q at the input of 38 
are approximately equal to the required real and imaginary 
parts of y y in equation (20), and thus the deconvolution of the 

20 multipath channel is achieved by replacing y by the j th 
component of the vector D^D.+^-lD^. 

Delay Estimation 

When the initial delay uncertainty x is not small compared 
with A, it is explicitly taken into account by the estimator 38. 

25 Denoting the initial uncertainty by x e -x p , the equation (19) 
is modified by replacing x by (x e -x p ), yielding 


D(r)=h 0 g c (x e -x p )+h 1 g c (x e -x 1 -x p )+. . . +h N g c (x.-x Ar x p )+n(x). (25) 

30 In (25) x p represents the negative of the signed integer part 
of (x/A) and x e is the fractional part of (x/A) such that 

T p =T e — £ 0 A; |x e |<A/2 

for some signed integer ko and with D selected above. The 
35 equation (25) may be rewritten as 

D(r)=h ogk ^)+h 1 g^(x e -x 1 )+. . . +h N g^{x-x N )+n{x) gkQ { x e)= 
gc(^-k 0 A) (26) 

Assuming that |k 0 |<N for some integer N, then ko is the 
40 integer k that minimizes the following index 




(27) 


45 


min f° +1)A \\BV-PMt4dT 
-N<k<N J^_ (io+1)A N II 


where h* is the channel response obtained on the basis of 
replacing ko by k in (26), 

g k (*HgkM,gk(^-'*i), ■ ■ ■ ,&(^-Xv)] r (28) 

and g^(x e ) is approximated by g c (-kA). In the estimation of 
h* for any integer k, the discretization and estimation pro- 
cedure of (20)-(22) is followed. Thus the estimator 38 
minimizes the following index 


60 


min 

-N<k<N 


M2 


Y 



(29) 


where g 0 * is the vector of dimension (q 1 +q 2 +l) whose 1th 
component is equal to g c (-(l+k)A), g k is equal to the j time 
shifted version of g*, and h* is the solution of equation (21) 
65 with g. replaced g k . 

If ko is the solution of the minimization in (29), then the 
overall channel estimate is given by h^lr 0 . Substitution of h 
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in (22) provides the coefficients of the inverse filter and the 
solution of (24) yields x, signal 46 in FIG. 1. Note that if it 
is known that the multipath error t is smaller than A/2 in 
magnitude, then this additional procedure is not necessary. 
However for k o >0, considerable error can otherwise ensue in 
the estimate of h. 

Multipath Phase Estimation 

With x denoting the estimate of the multipath delay error, 
one now solves equation (20) for h but with g- equal to 
g c (j A-x) by least squares estimation procedure. Denoting the 
resulting least squares solution by h^, then the multipath 
phase error estimate is given by the argument of the zeroth 
component of the estimated impulse response, i.e., 0=+arg 
(£</) which corresponds to signal 39 in FIG. 1. This error is 
compensated from the carrier phase-lock loop phase esti- 
mate in the following estimation cycle. 

Output 39 of estimator 38 is input to a loop filter 40, 
which may be a digital filter similar to the ones used in 
phase-locked loops. In its simple form it is a multiplier by 
a gain (1/G C T) where is G c the sensitivity of the numerically 
controlled oscillator (NCO) 42. Output 41 of the carrier loop 
filter 40 controls the phase of the carrier signal 43 (C-(t)) at 
the output of NCO 42. Signal 43 is inputted into phase 
shifter 44 which generates a jt/2 phase shifted signal 45 
(C^t)). Output 46 of estimator 38 is input to a code loop 
filter 47, which may be a digital filter similar to the ones used 
in delay lock loops or phase lock loops. In its simple form 
it is a multiplier by a gain (1/G^T) where G^ is the sensitivity 
of the NCO 49. Output 48 of the code loop filter 47 controls 
the phase (delay) of the output 50 of the code NCO 49. The 
output of code NCO is inputted to a code generator 51 whose 
output is the reference code a(t-x+Q 1 A) which is the binary 
pseudo random code and is the delayed version of the code 
a(t) appearing in the received signal 3 (S m (t). 

The delay t of this code is controlled by signal 48 at the 
input of code NCO 49. The code delay signal 52 is inputted 
to a tapped delay line 59 implemented by a shift register (not 
shown). The (Qi+C^+l) outputs of this tapped delay exem- 
plified by 52,53, and 54 are input to various correlators such 
as 8, 12, and 16, and to the correlators 23, 27, and 31. 

In the close-loop estimator, t is the estimate of the delay 
error x between the code phase of the direct component of 
the incoming signal 5 and the phase of locally generated 
signal 53 at one sampling time instance earlier. This error is 
corrected via the signal 48 at the input to the code NCO 49. 
With the recursive procedure in the limit the error t will 
approach zero and the code signal 53 will achieve lock with 
the direct path component of the incoming signal, thus 
eliminating the effect of multipath signal components 
present in signal 5. 

Similarly, 0 is the estimate of the phase error between the 
phase of the direct path component of the incoming signal 5 
and the local reference signal 43 (C/t)), and by generating 
signal 41 via lock loop filter 40, the control of the phase of 
carrier NCO 42 is such as to drive this error to zero. With the 
recursive procedure the carrier NCO output signal 43 will 
achieve phase lock with the direct path component of the 
incoming signal, thus eliminating the effect of the multipath 
signal components on the carrier phase tracking. The accu- 
mulated value of signal 41 multiplied by G C T or the accu- 
mulated value of signal 39 when the carrier loop filter 40 is 
a multiplier with gain (1/G C T), is the estimate of the differ- 
ence between the phase of the incoming signal direct path 
component and the initial carrier NCO signal phase. The 
value of signal 46 when the code loop filter is a multiplier 
with gain (1/G^T), is the estimate of the delay between the 
code phase of the incoming signal direct path component 
and the initial code phase of the code NCO output 50. 
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For the case of GPS signals, this delay difference multi- 
plied by the velocity of light is equal to the multipath error 
free estimate of the code phase range. Similarly the carrier 
phase difference multiplied by (k/ 2jt) where L is the carrier 
5 signal wavelength, is equal to the multipath error free 
estimate of the carrier phase range, also known as the 
accumulated delta range. 

SIMULATION RESULTS 

10 In this section some simulation examples are presented 
depicting the performance of the proposed multipath can- 
cellation algorithm. First some additional notations are 
introduced to present these results. Let T and T c denote the 
sampling period for the signal processing system and the 
15 code chip period respectively. Also let x d denote the nor- 
malized delay x d /T c and g be the sample signal-to-noise 
power ratio given by y=(P c T/N 0 ) where P c =A c 2 /2 is the 
signal power received by the direct or line-of-sight path and 
N 0 is the one-sided power spectral density of the receiver 
20 thermal noise. Further denote by a and 0 m the vectors 
consisting of the amplitudes and phases of the discrete 
channel response h in equation (26). One motivation behind 
this invention has been the GPS application to precision 
attitude determination of LEO (low earth orbit) spacecrafts. 
25 Therefore the achievable carrier phase estimation errors are 
interpreted in terms of the equivalent attitude pointing 
errors. It may be easily seen that for 1 meter antenna baseline 
the conversion factor between the carrier phase error in 
radians and the corresponding pointing error in arcmin is 
30 approximately equal to 109.4. Several simulation examples 
are presented below in terms of the notations introduced in 
this section. Recall that (q 1 +q 2 ) is equal to the number of 
discrete multipaths , (K 1 +K 2 +l) is the number of filter taps 
and g is the sample SNR. 

Example 1 

q 2 =0; q 2 =8; K a =12; K 2 =4; y=10 3 ;x^=0.5 
a=[l 0.2 0.5 0.8 0.9 0.7 1 0.2 0.9] 

40 0 m =[8 0 -0.5 0.7 1.2 0.8 -0.6 1.3 0.2 -1.1] 

FIG. 3 plots the ideal discriminator response g c (x), dis- 
torted response D(t) of (19), and the discrete version of the 
equalized response D eq (x). Clearly while the zero-crossing 
of D(t) is about 0.4 T c , the zero-crossings of both g c (x) and 
45 D eq (x) are equal to zero. Reducing the tap spacing to x^=0.1 
does not make any significant difference in the zero-crossing 
of the function D(t) as shown in FIG. 4, which plots the 
various discriminator functions for the case of x^=0.1. FIG. 
5 shows the convergence of the delay error as a function of 
50 the number of samples processed with an initial normalized 
delay of -0.498 chips. 

It can be seen that the steady state error without multipath 
correction will remain approximately equal to 0.4 T c while 
with the correction algorithm the error is approximately 
55 equal to 0.001 T c . Note also that with a 50 dB-Hz CNR at 
the receiver input this period of 100 iterations corresponds 
to only 1 sec of real time as per the definition of g. FIG. 6 
shows the carrier phase error expressed in arcmin for the 
GPS attitude determination application. Note that the initial 
60 error of more than 30 arcmin is reduced to about 0.1 arcmin. 

Example 2 

q 2 =0; q 2 =25; K 2 =35; K 2 =4; y=10 4 ; x^=0.1 
a=[l 0.8 0.5 0.3 0.9 0.2 0.5 0.3 0.7 0.1 0.05 0 0.08 0.1 0.03 
65 0.02 0.1 0.07 0.05 0.04 0.1 0 0.05 0.08 0.01] 

0 m =[O 0.5 0.7 -0.7 1 0.3 0.5 -0.6 1.5 1 0.8 0.2 0.1 0 1 -0.5 
-0.3 0.6 0.1 0.5 -0.8 -0.5 0.9 1.2 -0.6 1] 
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FIGS. 7 and 8 depict respectively the convergence of 
carrier phase and the code delay errors versus the number of 
samples. In this example the sample SNR is equal to 40 dB 
and thus corresponds to an update period of 0.1 sec for a 50 
dB-Hz CNR at the GPS receiver input. Note that the 5 
multipaths with delay greater than 9 D are weak in that their 
relative amplitudes are less than 0.1 of the direct path which 
makes the identification of these multipaths more difficult. 

As may be inferred from these figures a residual pointing 
error of less than 1 arcmin and a normalized delay error of 

0.005 chips is obtained in less than 50 iterations. Since the 
amplitudes of multipaths for delays exceeding 9 D are 
relatively small, one may consider ignoring these and thus 
select a filter order that is much smaller than q. FIG. 9 shows 
the resulting error performance when the filter order K a is 
selected to be 10. As is apparent from the figure the residual 15 
error, while much smaller compared to the initial error of 
about 30 arcmin, is much higher than that achieved with 
K a =35 filter case. 

Example 3 

F _ 20 

q a 0; q 2 =35; K-^35; 1^=4; y=10 3 ; x^=0.1 
a=[l 0.8 0.5 0.3 0.9 0 0.2 0.5 0.3 0.7 0.1 0.5 0.6 0.2 0 1 0.7 

0.5 0.2 0.3 0.5 0.7 1 0.4 0.8 0.3 0.7 1 0.7 0.8 0.6 1 0.5 0.3 

0 0.7] 

6 m =[0 0.5 0.7 -0.7 1 0.3 0.5 -0.6 1.5 1 0.8 0.2 0.1 0 1 -0.5 25 

-0.3 0.6 0.1 0.5 -0.8 -0.5 0.9 1.2 -0.6 1 0.7 -1 -0.7 0.8 

0. 6 0.5 -0.8 1 0.4 -0.9] 

FIGS. 10 and 11 plot the residual carrier phase and code 
delay errors respectively, for [g]=30 dB. As may be inferred 
from these figures, for this case of very severe multipath 
propagation, the pointing error is reduced from about 32 
arcmin to only 4.5 arcmin while the residual delay error is 
about 0.01 chips. 

It will be understood that the embodiment described 
herein are merely exemplary and that a person skilled in the 
art may make many variations and modifications without 35 
departing from the spirit and scope of the invention. For 
example, one such modification comprises sequential time 
sharing among a relatively small number of correlators to 
obtain correlation signals over a much larger range of delay 
offsets. This is achieved by sequentially offsetting the code 40 
NCO phase over a range of delay offsets and computing 
correlation signals for each such offsets instead of comput- 
ing all the correlation signals in parallel as described above 
in the preferred embodiment. Such sequential method may 
be attractive if the phase and delay of the input signal vary 45 
relatively slowly with time. All such variations and modi- 
fications are intended to be included within the scope of the 
invention as defined in the appended claims. 

We claim: 

1. A method for evaluating and compensating for multi- _ 
path signal error in code delay and carrier phase comprising 
the steps of: 

(a) receiving an input signal comprising a direct path 
signal and a plurality of multipath signals; 

(b) generating a first set of correlation signals by corre- 55 
lating the input signal with a known reference carrier 
signal and a plurality of reference code signals; 

(c) generating a real component of a discriminator func- 

tion by subtracting pairs of said first set of correlation 
signals; 60 

(d) generating a quadrature phase reference carrier signal 
obtained by shifting the known reference carrier signal 
by ninety degrees 

(e) generating a second set of correlation signals by 
correlating the input signal with the quadrature phase 65 
reference carrier signal and the plurality of reference 
code signals; 


(f) generating an imaginary component of the discrimi- 
nator function by subtracting pairs of said second set of 
correlation signals; 

(g) computing a carrier phase difference estimate between 
a phase of the known reference carrier signal and a 
carrier phase of the direct path component of the input 
signal, and also computing a code delay offset estimate 
between a selected reference code signal phase and a 
code phase of the direct path component of the input 
signal; 

(h) processing the carrier phase difference estimate using 
a carrier loop filter, and updating the phase of the 
known reference carrier signal to minimize the carrier 
phase difference; 

(i) processing the code delay offset using a code loop 
filter, and updating the selected reference code signal to 
minimize the code delay offset; 

(j) repeating steps (b) through (i) in real time such that 
said carrier phase difference and said code delay offset 
are maintained at values less than predetermined limits; 
and 

(k) outputting a signal corresponding to an estimated 
direct path component of said input signal. 

2. The method as recited in claim 1 further including the 
step of converting the received input signal from a carrier 
frequency to an intermediate frequency prior to generating 
said first set of correlation signals. 

3. The method as recited in claim 1 wherein the known 
reference signal is generated by a numerically controlled 
oscillator. 

4 . The method as recited in claim 1 wherein the plurality 
of reference code signals are generated by a numerically 
controlled oscillator, a code generator, and a tapped delay 
line. 

5 . The method as recited in claim 1 in which step (g) 
further comprises the steps of: 

(gl) generating a number of sampled multipath channel 
impulse response estimates corresponding to a number 
of uniformly spaced values of delay offsets using a least 
squares estimation algorithm; 

(g2) selecting one of the channel impulse response esti- 
mates based on minimizing the magnitude of a square 
of the difference between a complex valued sample of 
a discriminator function, whose real and imaginary 
components are respectively equal to said values of 
said first and second set of discriminator functions and 
sample values generated with the delay offset among 
said uniformly spaced values of delay offset and a 
corresponding channel impulse response estimate; 

(g3) generating a vector comprising the coefficients of an 
optimum deconvolution filter with said channel 
impulse response estimate and its shifted versions; 

(g4) convolving said complex valued sampled of said 
discriminator function with said coefficients of the 
deconvolution filter to generate a sampled version of an 
equalized discriminator function; 

(g5) estimating the code delay offset as the value for 
which the real part of the equalized discriminator 
function is equal to zero; 

(g6) generating a fixed multipath channel impulse 
response estimate based on said code delay offset 
estimate; and 

(g7) estimating the carrier phase difference as an argu- 
ment of the complex valued sample with an index equal 
to zero of said final multipath channel impulse response 
estimate. 
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6. A signal receiving apparatus which compensates for 
multipath signal error comprising: 

means for receiving an input signal comprising a direct 
path signal and a plurality of multipath signals, said 
direct path signal and said plurality of multipath signals 
each comprising a carrier component and modulated by 
a code component; 

means for processing said input signal in first and second 
sets of parallel correlators to generate first and second 
sets or correlation signals, said correlators comprising 
means for correlating said input signal with a reference 
carrier signal in the first set of correlators and with a 
ninety degree phase shifted signal with respect to said 
reference carrier signal in the second set of correlators, 
wherein each set of correlators correlate using delayed 
signals derived from a common reference code signal, 
said first and second set of correlators comprising at 
least three correlators; 

means for combining the correlation signals generated 
within each of said first and second set of correlators to 
generate values of a first and second discriminator 
functions, said values of the first and second discrimi- 
nator function comprising real and imaginary compo- 
nents of a complex-valued discriminator function; 
means for generating estimates of a phase error between 
a phase of the reference carrier signal and a carrier 
phase of the direct path component of the input signal; 
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means for generating estimates of a delay offset between 
a selected reference code signal phase and a code phase 
of the direct path component of said input signal; 
carrier loop filter means for compensating the phase of the 
reference carrier signal based on the estimate of the 
phase error to minimize said phase error; 
code loop filter means for compensating the phase of the 
reference code signal based on the estimate of the delay 
io offset to minimize said delay offset; and 

means for generating the direct path signal component of 
said input signal from said phase error and said delay 
offset. 

7. The signal receiving apparatus as recited in claim 6 
further comprising: 

a carrier numerically controlled oscillator, and 
a code numerically controlled oscillator, 
wherein said carrier loop filter means compensates said 
20 carrier numerically controlled oscillator to recursively con- 
verge the reference carrier phase with the direct path carrier 
phase, and said code loop filter means compensates said 
code numerically controlled oscillator to recursively con- 
25 verge the selected reference code signal phase with the direct 
path code phase. 



